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A recently obtained exact analytic solution to the wave-packet version of the pion 
- laser model is presented. In the rare gas limit, a thermal behaviour is found 
while the dense gas limiting case corresponds to Bose-Einstein condensation of 
wave-packets to the wave-packet mode with the smallest possible energy. 



1 Introduction 

The study of the statistical properties of quantum systems has a long history 
with important recent developments. In high energy physics quantum statis- 
tical correlations are studied in order to infer the space-time dimensions of 
the intermediate state formed in elementary particle reactions. In high energy 
heavy ion collisions hundreds of bosons are created. The correct theoretical 
description of their correlations is difficult. In this conference contribution we 
present a contradiction free treatement of multi-particle Bose-Einstein corre- 
lations for arbitrary large number of particles. 

2 Model Assumptions: 

A model system is described as 



P 



^Pnpn, (1) 



n=0 



the density matrix of the whole system is normalized to one. Here p n is the 

density matrix for events with fixed particle number n, which is normalized 

also to one. The probability for such an event is p n . 

The multiplicity distribution is described by the set {Pn} n=u - 

The density matrix of a system with a fixed number of boson wave packets 

has the form 

Pn = J da x ...doL n p n (ai, ...,a n ) ...,a n ){a u ...,a n \, (2) 

where \a\, ...,a n ) denote properly normalized n-particle wave-packet boson 
states. 
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The wave packet creation operator is given as 



at 



d 3 p 

(na 2 )'' 



»£i(p- 



7T i )+ia>(p)(t-t i ) . 

a T (p)- 



(3) 



The commutator 



(4) 



vanishes only in the case, when the wave packets do not overlap. 

Here on = / Ki,(Ji,ti) refers to the center in space, in momentum space, 
the width in momentum space and the production time, respectively. 

For simplicity we assume that on = (iTi,^i,a,to)- 

We call the attention to the fact that although one cannot attribute exactly 
defined values for space and momentum at the same time, one can define 
precisely the 7r,,£j parameters. 

The n boson states, normalized to unity, are given as 



ai, ... , a n ) = 



1 



(n) i=l 



4lo). 



(5) 



Here denotes the set of all the permutations of the indexes {1,2, ...,n} 
and the subscript <Ti denotes the index that replaces the index j in a given 
permutation from a^ n \ The normalization factor contains a sum of n\ term. 
These terms contain n different on parameters. Thus the further calculation 
with these normalized states seems to be extremely difficult. 



3 A New Type of Density Matrix 

There is one special density matrix, however, for which one can overcome 
the difficulty, related to the non- vanishing overlap of many hundreds of wave- 
packets, even in an explicit analytical manner. This density matrix is the 
product uncorrelated single particle matrices multiplied with a correlation fac- 
tor, related to the induced emission: 

p n (a u -,a n ) = -j^j ^flpiM^ (^2 LI ( ak \ a ^^j ■ ( 6 ) 



Normalization to one yields W(n). The weight factors describe induced emis- 



sion: 



p n {ai,...,a n ) (7 (n)fe=i . 

(overlap). (7) 



This is maximal if all are emitted with same ct\ 
p n (oti, ...,ai) n! 



[^i(ai)]" ~ W(n) 
and minimal if the overlap is negligible, 
p n (a 1 ,...,a„) 1 



(full overlap), (8) 



(no overlap). (9) 



11 Pi(«j) 

Thus we have wildly fluctuating weight. E.g. for 800 pions the induced emis- 
sion weight fluctuates between [1, 800!] ~ [1, 10 1977 ]. With this special density 
matrix one can proceed with the calculations. 

3.1 Single- Particle Density Matrix : 

For the sake of simplicity we assume a factorizable Gaussian form for the 
distribution of the parameters of the single-particle states: 

Pi{a) = p x (0 p p (tt) 5(t - t ), 

PAO = — ^— r cM-e/(2R 2 )), 
(2t:R z )2 

PpW = - ^ 3 exp(-7r 2 /(2mT)). (10) 
(zirml ) 2 

These expressions are given in the frame where we have a non-expanding static 
source at rest. 

A multiplicity distribution when Bose- Einstein effects are switched off (de- 
noted by Pn^), is a FREE choice in the model. We assume independent emis- 
sion, 

= ^exp(-no). (11) 
n! 

This completes the specification of the model. 
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4 Solution of Recurrences: 



The probability of finding events with multiplicity n, as well as the single- 
particle and the two-particle momentum distribution in such events is given 
as 



Pn 



2-1 



N. 



(»), 



ki,k a ) = £)Z) ^[G ro (k 1 ,k 1 )G i _ m (k 2 ,k 2 



(12) 



G m (k 1 ,k 2 )G i _ m (k 2) k 1 )] 



(13) 



In rem recurrence relations were given for the construction of the func- 
tions G m (ki,k 2 ). In refstl, □ an explicit analytic form was obtained for these 
functions. We present now this method. 

To arrive to this solution we introduce the following auxiliary quantities: 



7± = - (1 + a ± VI + 2x) 



K p a T . 



Using the notation 

or = 2mT p , 



T 



2m' 



Rz 



R A 



mT 

(7 (Jrp 



(14) 



(15) 



we arrive to the formulae of the plane- wave pion laser model of Pratt, ref.u, if 
we replace R in ref. u with R p and T in ref. u with T p . The general analytical 

(16) 



solution of the model is given through the generating function of p n , 

oo / oc \ 

G(z) = ]Tp„z" = exp £ C„(z" - 1) , 

n=0 \n=l J 

where C n is the combinant 01 of order n, 



n 



n 



7+ - 7- 

and the general analytic solution for the functions G„(k 1; k 2 ) reads as: 



G„(ki,k 2 ) = j„c 



ki-7jk ; 



1 7+ - 7- 
4 1+ - 7- ' 



(17) 

(18) 
(19) 
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The detailed proof that the analytic solution to the multi-particle pa/aye-packet 
model is indeed given by the above equations is described in refs.El,l3. 

With our analytic solution we decreased the algorithmic complexity of the 
problem of n symmetrized boson system from the original n\ to 1. 

5 Multiplicity Distribution: 

The mean and the second factorial moments are defined as follows 



(n) = 5>p„ = Y^iCi, (20) 

n i—1 

oo 

i(n-l)) = ^n(n-l)p„ - (n) 2 + ^ * (i - 1) Q. (21) 



i=2 



The large n behavior of the multiplicity distribution, p n , depends on no/7?. 
This multiplicity distribution has an interesting property. To see this, we define 
the quantity n c as 



7+ 



1 + x + \/\ + 2a 



(22) 



One can show that (n) is finite, if no < n c and (n) is infinite, if no > n c Hence 
n c is a critical value for the parameter no- 



6 Solution for the Inclusive Distributions: 

In the previous sections we obtained results the multiplicity distribution and 
for exclusive momentum distributions. To obtain the inclusive distribution we 
introduce the auxiliary quantity 

00 

G(k l5 k 2 ) = ^Gi(ki,k 2 ). (23) 

i=i 

( Higher order Bose-Einstein symmetrization effects are negligible, if the first 
term dominates the above infinite sum, i.e. if G(ki,k2) = Gi(ki,k2).) 
Now the averaging over the exclusive distributions can be performed, 

00 

7V x (k) = £> n JV{">(k) = G(k,k), (24) 

n=l 
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and the two-particle inclusive correlation functions can be evaluated as 

iV 2 (ki,k 2 ) 



G 2 (k l5 k 2 ) = 



JVi(ki)JVi(k 2 ) 

G(ki,k 2 )G(k 2 ,ki) 



1 + 



G(k 1 ,k 1 )G(k 2 ,k 2 



(25) 



This is an exact result, obtained without the random phase approximation and 
without assuming that the number of sources is large. However, this result is 
valid only when Bose-Einstein condensation does not play a role, i.e. when 

n < n c = 7 + ' . 

7 Rare gas limiting case: 

For large source sizes or large effective temperatures, x » 1 we have : 

|(n-l) 



G„(ki,k 2 ) oc — 

x 



G„ 



exp 

|(l-n) 



2<Jrp 



kD-^ck.-k,) 2 



, (26) 



(27) 



From this equation we can see that the effective temperatures and the effective 
radii are decreased by a factor of f/n for n-th order symmetrization effects. 
The multiplicity distribution is shifted to high multiplicities: 



Pn 



— - exp(-n ) 
n! 



f + 



n(n — f ) — tjq 
2(2x)2 



(28) 



The momentum distribution is obtained in this approximation by an ex- 
pansion in the small parameter e = (2/x) 3 ^ 2 . In this approximation we get 



JVi(k) 
A^k) 



n 



(ttct t ) 
n 



2 s a ex P 



(ttct t ) 



2 ^ 2 eX P 




The single-particle inclusive and exclusive momentum distributions are 
enhanced at low momentum. The influence of the wave packet size on the pion 
multiplicity distribution is shown in Fig.l. The effect of the wave packet width 
on the momentum distribution is displayed on Fig. 2. 
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Figure 1: As the size of the wave packets is changed, the overlap changes and stimulated 
emission results in larger multiplicities. 
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8 Correlation functions: 

In the highly condensed x << 1 and n >> n c Bose gas limit a kind of 
lasing behaviour_and an optically coherent behaviour is obtained, which is 
characterized byl3 the disappearance of the bump in the correlation function: 

C 2 (k 1 ,k 2 ) = l. 

On the other hand, in the rare gas limit we get 

Ca(ki,ka) = 1 + exp (-i? 2 Ak 2 ) . (31) 



.1 Exclusive Correlation Functions: 



In a highly condensed limiting case, x << 1 and tiq >> n c , the exclusive and 
inclusive correlation functions coincide, C*2(ki,k2) = 1. 

In the rare gas limit, x » 1, the exclusive and inclusive correlations are 
different, and the exclusive correlation has the form 



^ 2 (n) (ki,k 2 ) 



n(n-l) AT 1 (n) (ki)iV}™)(k2) 
1 + A K exp (-i? K „ Ak 2 - i? 2 Ak 2 



(32) 



k 2 ), Ak, = Ak - K(Ak-K)/(K • K) and Ak Q 



where K — 0.5(ki 
K(Ak-K)/(K-K), 

The momentum dependent parameters are given as follows: 



A 



K 



R 2 



K. 



1 



Rl 



1 



1 



2( 5 / 2 ) exp 



R 



K,o 



R 



(2x)i 
n 

J'2 



R 



K 2 

— exp 

0~rp 



\Pi exp 
K 2 



K 



K 2 



2 \ 1 



(n + 2)R. 



„2 
(J T 



(33) 
(34) 
(35) 



The dependence of the parameters of the correlation functions on the mean 
momentum of the two pions, K, is shown on Fig. 2. The influence of the 
wave packet size on the two particle correlation is-illustrated on Fig. 3. (These 
limiting cases are discussed in more details in refill.) 
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9 Highlights: 

In this paper we presented a consequent quantum mechanical description of 
the correlations caused by the multi-particle Bosc-Einstcin symmetrization of 
a system of large number of bosons. 

We introduced the overlapping multi-particle wave-packet density matrix 
describing stimulated emission. We reduced the algorithmic complexity of 
the description of the n boson states from n! to 1 by obtaining an cxplicitc 
analytical solution for the problem. 

We have shown, that the radius and intercept parameters depend on the 
mean momentum of the two pions even for static sources, due to multi-particle 
symmetrization effects. 

We have found an enhancement of the wave-packets in the low momentum 
modes, due to multi-particle Bosc-Einstein symmetrization. When all pions are 
in the lowest momentum state, the system may be interpreted as a pion-laser. 
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Multi-Particle Symmetrization Effects 
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Figure 4: Momentum-dependent reduction of the intercept parameter Ak; the side- wards 
and the outwards radius parameters, i?K,s an <i Rk.,o from their static values of 1 and R p , 
respectively. 
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